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Abstract 

We analytically study the linear propagation of arbitrarily shaped light-pulses through an ab- 



^ , sorbing medium with a narrow transparency-window or through a resonant amplifying medium. 



We point out that, under certain general conditions, the pulse acquires a nearly Gaussian shape, 
irrespective of its initial shape and of the spectral profile of the line. We explicitly derive in this 
case the pulse parameters, including its skewness responsible for a deviation of the delay of the 
pulse maximum from the group delay. We illustrate our general results by analysing the slow-light 
experiments having demonstrated the largest fractional pulse-delays. 
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I. INTRODUCTION 



The principle underlying most of the slow-light experiments is to exploit the steep normal 
dispersion of the refractive index associated with a pronounced peak in the transmission of 
the medium and the correlative reduction of the group velocity. The situation where the 
resulting time-delay of the light-pulse is large compared to its duration and can be controlled 
by an external laser field is of special importance for potential applications, especially in 



the domain of high-speed all-optical signal-processing. Harris and his co-workers 
opened the way to such experiments by exploiting the phenomenon of electromagnetically 
induced transparency (EIT) allowing one to create a narrow transparency-window in an 
otherwise optically thick atomic vapour. Using a true-shape detection of the pulses, they 
demonstrated propagation velocities as slow as c/165 and group delays Tg as long as 4.6 Tin 
where c and rj„ are respectively the velocity of the light in vacuum and the full width at 
half-maximum (FWHM) of the intensity-profile of the incident pulse. Much slower velocities 
have been attained in subsequent EIT experiments (for reviews see, e.g., and 
in experiments involving coherent population oscillations |6] or other processes to induce 
a transparency-window in an absorbing medium. It is however worth noticing that only 
few of these experiments, all using EIT, have succeeded in giving direct demonstrations of 
fractional delays Tg/rin exceeding unity . Theoretical discussions on the maximum 

time-delays attainable in such experiments can be found in ^ , y, Q| . A different way to 
achieve a system with a controllable transmission-peak is to opticallv induce a resonant 
gain .n a n,ediu„> Q. MtiaUv — d by Gauthie. Q. the a^nge.ent 

involving stimulated Brillouin scattering [13|, [IJ, [151, [161] seems particularly attractive from 
the viewpoint of the above mentioned applications. The Brillouin gain is indeed directly 
implemented on an optical fibre and there are no severe constraints in the choice of the 



operating wavelength. The group delay Tg has already been controlled on a range of 3.6 rj„ 
by this technique Note that preliminary experiments using a Raman fibre amplifier 

have also been achieved [l7|. 

The purpose of our paper is to provide analytical results on the propagation of arbitrarily 
shaped pulses, the central frequency of which coincides with that of a pronounced maximum 
in the medium-transmission. We examine more specifically the case where the resulting time- 
delays of the pulses are large compared to their duration. Our study applies in particular but 
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not exclusively to the above mentioned systems. Our approach follows in part and extends 
that of Bukhman with a special attention paid to the connection of the theoretical 
results with the experiments. 



II. GENERAL ANALYSIS 

We denote by ej„(t) and Coutif) the slowly- varying envelopes of the incident and trans- 
mitted pulses and by EiniVt) = emit) exp{—iQt)dt and EoutiS^) their Fourier transforms. 
The slow-light medium is characterised by its impulse response h{t) or by its transfer func- 
tion H{Q), Fourier transform of h{t). The input/output relation or transfer equation reads 
en„.t (t) = h{t)®ein{t) in the time-domain or Eout{^) = H{Q)Ein{Q) in the frequency-domain 
(l9l |. We assume that the incident pulse has a finite energy, that it is not chirped (ej„(t) real 
and positive) and that h{t) is also real. The local response of the medium is characterised 
by the complex gain-factor T{Q) = In whose real part F{Q) and imaginary part 

are respectively the logarithm of the medium amplitude-gain and the induced phase 

shift. The condition imposed to h{t) implies that H{—Q) = H*{Q) and thus that F{Q) and 
$($1) are respectively even and odd functions of Q. This has the advantage of eliminating 
the lowest-order pulse-distortions resulting from the gain-slope and from the group velocity 
dispersion at the frequency c^o of the optical carrier (fi = 0). Moreover the medium is then 
entirely characterised by the single real function h{t). In order to have simple expressions 
we use for emit) a time origin located at the pulse centre-of-gravity and for eout{t) a time 
origin retarded by the transit time at the group velocity outside the frequency-domain of 
high-dispersion (local time picture). The time delays considered hereafter are thus only 
those originating in the high-dispersion region. 

General properties of the transmitted pulse can be derived by Fourier analysis. Let x{t) 
be any of the real functions em(t), h{t) or eout{t) and X{Q) its Fourier transform. We 
remark that X{0) = x{t)dt and, following an usual procedure in probability theory |2ol |. 
we characterise X{Q) by its cumulants k„, such that 



For H{Q), we see that the cumulants are simply related to the coefficients of the series 
expansion of T{Q) in powers of —iQ and, in particular, that Hi coincides with the group 




(1) 
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delay Tg = \n=Q- We incidentally recall that, due to the causality principle, Tg can be 
related to the gain profile 
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22l |. Within our assumptions, this relation reads 



where Hq = H{0) is the amplitude-gain of the medium at the frequency of the optical carrier. 
This confirms that large group delays are achieved when the gain has a pronounced 

maximum at f2 = 0. We have then k,2 > 0. 

Coming back to the general problem, we characterise the time function x{t) by its 
area S = x{t)dt and its three lowest order moments, namely the mean value (t) = 
-| x{t)dt, the variance = /-^(^ ^ {t))'^x(t)dt and the 3'^'^ order centred moment 
a = ^ J^(t— {t))^x{t)dt . We recognise in {t) (resp. a) the location of the centre-of-gravity 
(resp. the rms duration) of the function x{t). Its asymmetrvmay be characterised by the 
dimensionless parameter ^ = a/a^, the so-called skewness j20|. For a Gaussian function, 
cr = r/ (2Vln2) where r is the FWHM of the energy profile x'^{t). An important result [2^ 
is that the moments (t) , and a of x{t) are equal to the cumulants Hi, k,2 and of 
Moreover the transfer equation immediately leads to the relations Eout{^) = -f^o-£'m(0) and 
i^n,out = i^n,m + i^n wherc, as in all our paper, the indexes in, out and the absence of index 
respectively refer to the incident pulse, the transmitted pulse and the transfer-function or 
impulse response of the medium. With our choice of time origin, (tin) = 0. By combin- 
ing the previous results, we finally obtain the four equations Sout = H^Sin, {tout) = Tg, 
^out — ^in + o"^ > o'in '^out = o-in + 0,. In the studies of the linear pulse propagation, the 



out 

first equation which relates the areas of the transmitted and incident pulses is known as the 
13]. The second equati 

^ y. 



area theorem |22|. The second equation expresses that the time-delay of the pulse centre- 
of-gravity equals the group delay |2^. The two last ones specify how the rms duration and 



the asymmetry of the incident pulse are modified 
provided that the involved moments are finite jl9|. 



yy the medium. All these results are valid 



III. ANALYTIC EXPRESSIONS OF THE MEDIUM IMPULSE-RESPONSE 

In order to obtain a complete information on the shape and the amplitude of the trans- 
mitted pulse, we obviously have to specify the complex gain-factor r(fi) of the medium. For 
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the medium with a resonant gain, it reads 



V{SI)=pn{^)G/2- A/2 (3) 

where Pn{^) is the normahsed complex profile of the gain-line (pAr(O) = 1), G is the gain 
parameter for the intensity ^| and A stands for the attenuation introduced to reduce the 
effects of the amplified spontaneous emission and/or to normalise the overall gain of the 
system. G = QqL where goi^resp-L) is the resonance gain-coefficient (resp. the thickness) of 
the medium. The previous expression of also holds for an absorbing medium with a 
transparency-window when the absorption background is assumed to be infinitely wide. We 
then get T{Q) = —[1 — f pj^{Q)]aoL/2 where ao is the background absorption-coefficient, 
/ < 1 specifies the depth of the transparency- window [9] and Pn{^) is the normalised 
complex profile of the line associated with the transparency- window. By putting G = f a^L 
and A = a^L , we actually retrieve the expression of obtained for a gain medium. In 
both types of experiments, G and A are generally comparable in order that the resonance 
gain Hq = e*^/^""^/^ is close to 1 or, at least, does not differ too strongly from 1. Anyway the 
intensity-transmission on resonance exceeds its value far from resonance by the factor e*^. 

To go beyond, it seems necessary to explicit the profile pn{^)- We first consider the 
reference case where Pn{^) is associated with a Lorentzian line. It then reads Pn{^) = 
1 /{l+iQ/'-f) where 7 is the half- width of the line and we immediately get k„ = Gn\ / (27") 
with in particular ki = Tg = G/{2'y), K2 = cr^ = G/'y'^ and thus Tg = a\/G/2. The last 
relation shows that achieving substantial fractional delays Tg/a requires that G ^ 1. A 
quite remarkable property of the Lorentzian case is that the impulse response has an exact 
analytical expression. This result has been obtained by Crisp in a general study on the 



propagation of small-area pulses in absorbing and amplifying media 



23[ but it can easily be 



retrieved from H{Q) by using standard procedures of Laplace transforms j^]. One get 



hit) = e-^/^5(t) + e-^/%G ^^^^^ e-*f/(t) (4) 

where S{t), Ii{u) and U{t) respectively designate the Dirac function, the order modified 
Bessel function and the unit step function. The 1'^* term hi{t) in h{t) results from the con- 
stant value e~^/^ of H{^1) far from resonance. This part of the response is instantaneous 
in our local time picture and only the 2"''^ term hd{t), directly associated with the trans- 
mission peak, contribute to the delay. The areas of hi{t) and hd{t) are respectively e""^/^ 



and Hq — e that is in a very small ratio (~ e "-^/^ ) for the large values of G required 
to achieve substantial delays (see above). The effect of the instantaneous response then be- 
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FIG. 1: Analytical form of the impulse response for a Lorentzian line. From the left to the right 
the gain parameter G (resp. the fractional delay Tg/a) is 4, 9, 16, 25, 36, 49 and 64 (resp. 1, 1.5, 
2, 2.5, 3, 3.5 and 4). The horizontal (resp. vertical) unit is a (resp. HQ/aV2TT ). 



comes negligible. Fig^ shows the delayed response obtained for increasing values of the gain 
and thus of the fractional group-delay. We see that the curves, first strongly asymmetric, 
become more and more symmetric as G increases and that the location of their maximum 
then approaches the group delay. They have a discontinuity H^e^'^/'^'yG /2 at t = 0, the rel- 
ative amplitude of which becomes negligible when G ^ 1. From the asymptotic behaviour 
of Ii{u) we then get 
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hAt) f 1 - HI exp (-4 ) (5) 



(TV27r V 4rgy V 
with 6 = t — Tg. The maximum of /^^(t) occurs at the instant Tg — At with At ^ 3cr^/(4rg). 
When G ^ oo. 

This Gaussian form is that of the normal distribution derived by means of the central limit 
theorem in probability theory. This theorem can also be used for an approximate evaluation 
of the convolution of n deterministic functions It applies to our case by splitting the 
medium in n cascaded sections, h{t) being then the convolution of the impulse responses 
of each section. According to this analysis, one may expect that the normal form h^'^\t) 



is universal. From the frequency viewpoint, it originates in the fact that, when G ^ 1 , 
the transmission peak is roughly a/G times narrower than the line. In the region where the 
relative gain / Hq is not negligible, the curves vs VL (phase-shift) and FiVL) vs 

f2 (line-profile) are well approximated respectively by a straight line and a parabola (FiglSj). 
This means that only the first two cumulants Ki = Tg and K2 = cr^ play a significant role. 
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FIG. 2: ^{^) and F{0,) associated with a Lorentzian line (full line) and a Gaussian line (dashed 
line). The physical parameters are chosen in order that Hq = 0.6 and Tg/a = 3 in both cases. The 
resulting gain profiles |-ff(r2)| are indistinguishable at the figure scale. The frequency unit is 1/a 
(angular frequency). 



We then get H{n) ^ H^'^\Q) = Hoexp^-iQTg-a'^n'^ /2) and thus h{t) ^ h^'^\t) irrespective 
of the line-profile. This confirms the universality of the normal form of h{t) when G — » 00. 
In fact, h^'^\t) is a good approximation of the exact impulse- response for the gains currently 
achieved in the experiments. Fig|21 shows the result obtained with Lorentzian and Gaussian 
line-profiles when Tg = 3a. In the second case, Pn{^) = exp(— ^7^/7^) — 2iD(yL/'~^) / ^/tx with 
D{u) = e~^^ J^e^^dv 24| and the first cumulants read ni = G/{'yy/n), K2 = and 
K3 = 4G / {'-f^ ^/^T) . The parameters A, G and 7 are chosen such that Hq, Tg and a have 
the same values in both cases. Though the line-profiles are quite different (see Fig|21), the 
impulse responses are both close to the normal form h^'^\t). Similar results (not shown) are 
obtained with other line-profiles, including the FIT profile (see hereafter). 

When the gain parameter is large but not very large, a better approximation of the 
impulse response is obtained by considering the effect of the 3'"'^ cumulant K3, equal to 
the asymmetry parameter a of h{t). Provided that this effect may be considered as a 
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small perturbation, H{Q) ^ H^^\Q) ^ {1 + ik^Q^ /3\) H^'^\fl). From the correspondence 
iQ ^ d/dt, we finally get 

h{t)^h^'\t)^(i-^y(^\t). (7) 

This result generalises that obtained in the Lorentzian case where a = 3G'/7^ and a/ (2cr^) = 
3/{4Tg). In the Gaussian case, we find a/(2(j^) = 2/{7CTg), a value not far from the previous 
one. This explains why the two impulse responses are very close (see Figl^I). Anyway they 
are very well approximated by h^^\t) in each case. Quite generally the maximum of h^^\t) 
occurs at Tg — At with At ^ a/(2cr^) = where C, is the skewness of h{t). In all the above 
calculations we have implicitly assumed that At is small compared to a. This implies that 
1^1 ^ 2 but we checked that h^^\t) keeps a fairly good approximation of h(t) for skewness 
up to 1. 
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FIG. 3: Comparison of the exact impulse-responses for the Lorentzian profile (full line) and the 
Gaussian profile (dashed line) to the normal form h^'^\t) (dotted line). The fit by the improved 
forms h^^\t) (not shown for clarity) is nearly perfect. Parameters as in FigJ^I Units as in Fig^ 

IV. NORMAL FORM OF THE TRANSMITTED PULSE 

The impulse response being known, the envelope of the transmitted pulse is given by 
the relation eout{t) = h{t) ^ ein{t) and will generally differ from that of the incident pulse. 
However the distortion will be negligible if the duration of h{t) is small compared to that of 
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the pulse. We then get h{t) ^ 5{t — Tg) h{t)dt and thus eout{t) ~ Hoein{t — Tg). Since we 
are interested in the situations where the time delay is large compared to the pulse duration, 
we obtain the double condition a <^ ^ Tg which can only be met with extremely large 
gain parameters. Taking for example a = am/ 7 and Tg = that is Tg = 4:9a, we get 
G ~ 9600 in the Lorentzian case. Fig. 4 shows the results obtained with these parameters. 
As expected the pulse distortion is small, even in the sensitive case of a square-shaped pulse. 
Note that the gain parameter considered is not unrealistic. It is comparable to that used by 
Harris and co-workers in their pioneering EIT experiment where G ^ A ^ 6000 ^. 
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FIG. 4: Propagation of a square-shaped (full line) and a Gaussian-shaped (dashed line) light-pulse 
with large delay and low distortion. The parameters are chosen in order that a = (Jin/7, Tg = 7ain 
and Hq = 1. The envelopes of the input pulses are given for reference. The time unit is the common 
rms duration o"j„ of the two incident pulses. As expected, the distortion of the Gaussian-shaped 
pulse is negligible and there is only a slight softening of the rise and of the fall of the square-shaped 
pulse. 



However most of the direct demonstrations of large fractional pulse delays have been 
achieved with smaller gain parameters, typically ranging from 10 to 100, and with incident 
pulses whose rms duration am is comparable to and often smaller than a. Substantial 
pulse-reshaping is then expected. Suppose first that the normal form /i*^^) {t) provides a good 
approximation of the impulse response and that the incident pulse is Gaussian-shaped. We 
have then em{t) = exp(— t^/2(T^) and Coutit), convolution of two Gaussian functions, is itself 
Gaussian. It reads Coutif) = Hoi^am/ aout) /'^'^lut) ^ where a^^^ = af^ + a"^ . The effect 

of the medium on the light-pulse is simply to delay its maximum exactly by the group delay, 
to broaden it by the factor \/ 1 + a/^ /af„ [1, 9] and to modify its amplitude accordingly in 
order to respect the area theorem j23|. Since this point is often overlooked, we stress that 
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ihe broadening mechanism radically differs from that occurring in standard optical fibres 
2^ • It originates in the 2"'^ order gain-dispersion instead of in the group- velocity dispersion 
and the pulse envelope keeps real (no phase modulation or frequency chirping). In fact, 
provided that crj„ be smaller than or comparable to a and that l^j^l < 1, Coutit) is well 
approximated by a Gaussian function whatever the shape of the incident pulse is. This is 
again a consequence of the central limit theorem, the response eoutit) being obtained by an 
extra convolution added to those used to build h(t). The conditions on am and \C,in\ originate 
in the requirement that all the terms to convolute should have moments of the same order of 

(2) (2) 

magnitude. We then obtain Coutit) ~ elutit) where elJtit) has the normal (Gaussian) form 

e?lit) = i/o— exp (--^) . (8) 

This result extends the previous one and shows that incident pulses having different shapes 
but the same area and the same variance 0"^ are reshaped in the medium to give 
approximately the same Gaussian-shaped pulse (FiglHl). From an experimental viewpoint, 
the dramatic reshaping of a square-shaped pulse has been clearly demonstrated (but not 
commented on) by Turukhin et al. [8] in their EIT experiment in a solid (see their figure 2c 
for probe detuning). Pulse reshaping is also apparent in the Brillouin scattering experiment 
by Song et al. [l^ where a fiat-topped pulse is actually transformed in a gaussian-like pulse 
(see their figure 4 and compare the shapes obtained for gains OdB and 30dB). 

A more precise approximation elJ^^it) of eout{t) can be obtained by taking into account 
the effect of the S*"*^ order cumulants. Using the approach already used to determine h^^\t), 
we get 

eS(t) ^(l-^) e^it) (9) 



with ttout = ttin + d- When the incident pulse is symmetric (aj„ = ) as in most experiments, 
the skewness ^out of the transmitted pulse reads C,out = a/ (o"^ -|- o-fnf^"^ and the pulse maxi- 
mum occurs at Tg — Atout with Atout = o./ [2 (a^ + crfj^)]. Since \C,out\ < \(,\, the transmitted 
pulse is closer to a normal form than the impulse response of the medium. The previous 
results hold without restriction to the value of (7^ when ej„(t) is Gaussian. In the case of a 
Lorentzian line-profile, we easily get 

3 

= 27(l + af„7VG)- 

We have compared the theoretical delay of the pulse maximum, namely Tg — Atout, with 
the delay actually observed by Okawachi et al. in their Brillouin scattering experiment 
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FIG. 5: Example of pulse reshaping and broadening. The square-shaped (full line) and Gaussian- 
shaped (dashed line) incident pulses originate nearly identical transmitted pulses, respectively close 
and very close to the normal form e^uti^) (dotted line). The time unit is cjj„ and the parameters 
are such that a = am, Tg/a = 4 and Hq = 1. With this choice of Hq and a, the transmitted pulses 
have the same area that the incident ones (area theorem) and a rms duration \/2 times larger. 

Ql . Figini shows this delay as a function of the gain parameter G for two values of the 
pulse duration, respectively rj„ = 63 ns (crj„ ^ 38 ns) and rj„ = 15 ns {(Xin ~ 9 ), with 
7 = 0.22 ns~^ (7 is the half of the full Brillouin linewidth Tb)- Without any adjustment 
of parameters, our analytical results satisfactorily fit the observations. Note that the shifts 
Atout are negligible for the longer pulse but significant for the shorter one. 



V. EFFECT OF THE TRANSMISSION BACKGROUND 

In the previous calculations, we have not taken into account the effect of the instantaneous 
part hi(t) = e~'^^'^5{t) of the impulse response, arguing that its area is small compared to 
that of the delayed part. As a matter of fact, hi{t) originates a contribution e~'^/^ej„(t) to 
^out{t), the amplitude of which is roughly Q^^'^Oinj Oo^t times smaller than that of the main 
part and is thus actually negligible in every case of substantial delay (G ^ 1, Oout and a in 
of the same order of magnitude). We should however remark that this result lies on the 
assumption that the transmission peak is put on an uniform background. 

We will now examine the case where the transmission background is not uniform. This 
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GAIN PARAMETER (G) 



FIG. 6: Comparison of the delays of the maximum of the transmitted pulse observed in a Brillouin 
scattering experiment [l^ (filled squares and circles) with our analytical predictions (full lines). 



happens in all the experiments where a transparency-window is induced in a absorption- 
profile of finite width, in particular in the EIT experiments. As an illustrative example we 
consider the simplest A arrangement with a resonant control field. From the results eriven 

n 

in 1^, we easily get 

where (resp. 7ca -C 76a) is the coherence relaxation-rate for the probe transition (resp. 
for the forbidden transition), VL^ is the modulus of the Rabi frequency associated with 
the control field and A = a^L ^ 1 is the resonance optical thickness in the absence of 
control field. The control field makes the resonance gain rise from e'^/"^ ^ to Hq = 
exp (— v47ba7ca/ [2 {jbalca + ^s/4)]) and a good transparency is induced when is larger 
than or comparable to 'jbalcaA. The width of the transparency-window (oc ils/ ) is then 
much smaller than that of the absorption background (cx 'jbaVA). 

Without any approximation, the partial fraction decomposition of T{Q) allows us to 
write the transfer-function of the medium as a product of simpler functions, namely 
H{Q) = Hi{Q)H2{^l) with Hj{Q) = exp[Cj/ {iQ + 'jj)]. According to the control power, 
the parameters Ci, C2, 71 and 72 are real or complex. When Qg < {iba — lea), 7i and 
72 are real and positive whereas Ci and C2 are also real but of opposite sign. The EIT 
medium is then equivalent to a medium with two Lorentzian lines both centred at f2 = 0, 
respectively an absorption-line and a narrower gain-line. It is also equivalent to a cascade 
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of a gain medium and an absorbing medium. When Qs > {jba — lea), all the parameters 
are complex with 72 = 7^ and C2 = CI . The two lines are now located at Q = ±Im(7i). 
They have the same intensity and the same width, but they are hybrid in the sense that, 
due to the complex nature of Ci and C2, their absorption and dispersion profiles are both 
the sum of an absorption-like and a dispersion-like profile. We incidentally note that the 
parameters used to obtain the figure 8 in ^ correspond to such a situation. In all cases, the 
impulse responses associated to Hi{Q) and H2{^1) have analytical expressions l^l and 
the impulse response h{t) of the medium is their convolution product. This general analysis 
is satisfactory from a formal viewpoint. It provides some physical insight into the EIT mech- 
anisms but is not really operational to determine the shape of the transmitted pulse. From 
this viewpoint, a fruitful approach consists in exploiting the fact that the medium is opaque 
except in the narrow region of induced transparency and in the far wings of the background 
absorption-line (Fig|7j). In the first region, H(^l) is well approximated by the forms H^'^\fl) 
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FIG. 7: Gain profile in an EIT experiment and expanded view of its central part (upper scale). In 
the wings as in the transparency window, the exact gain (full line) is scarcely distinguishable 

from the approximate form \H^'^\^}) + Hoff{0,)\ (dashed line). Remind that |//(i7)| ^ 1 in the far 
wings. The parameters are A = 63, 7ba/27r = 5 MHz, jca/jba = 2.2 10^^ and Q^s/lba = 0.73 . The 
frequency unit is 10^ Rd/s (angular frequency). 



or, if necessary, H^^\Vl) , obtained by keeping only the 2 or 3 first cumulants of HiVt) as 
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in the case of an uniform background. We only give here the simphfied expressions of these 
cumulants when <^ 'Jba and ^ jcalba (conditions of good induced transparency). We 
then get = Tg ^ 2A^ka/^l = ^ ^QA^L/^t and = a ^ 48^7^,(472^ - n2)/j^6 
with Hq ^ exp {—2A'ji,a'jca/^'i)- In the far wings \Q\ ^ Qg and H{Q) ^ Haff{Q), where 
Hoff{^l) = exp(— yl/[2(l +iQ/'yha)]) is the transfer-function when the control field is off. 
Finally we get the relation H{n) ^ H^p^Q) + Hoff{fl) with p = 2 or 3, valid at every 
frequency. Figi?! obtained for typical physical parameters, shows that \H'-'^\il) + Hoff{^)\ 
already provides a good approximation of the exact gain. Now reduced to a simple sum in- 
stead of a convolution product, the medium impulse-response reads h{t) = h^P\t) + hoff{t) 
where hoff{t), associated with a Lorentzian absorption-line, has an analytical expression 
I23I 1. As in the case of a gain-line, this expression can be retrieved from Hoff{^l) by using 
standard procedures of Laplace transforms [2^ . It reads 

hoffit) = 6{t) - 7..A '^^^^g2?^ exp (-7..t) U{t) (12) 

where Ji{u) designates the ordinary 1** order Bessel function. The envelope eout{t) of the 
transmitted pulse will be thus the sum of two terms. The first one is the approximate 
solution ef}^f-{t) obtained by the cumulants procedure. The second one reads eoff{t) = 
hoff{t) ®ein{t). It is worth remarking that hoff{t) is rapidly oscillating (characteristic time 
oc l/Ajha) and that its area is very small {f^^ hoff(t)dt = Hoff{0) = e"^/^ ). This entails 
that eoff{t) will have a negligible amplitude (oc e"^/^) when ej„(t) is smooth enough so 
that the far wings of its Fourier spectrum do not overlap those of the absorption-line. FiglHl 
shows the result obtained with a Gaussian-shaped incident-pulse. A , 7ba and crj„ being 
given, we have chosen the other parameters in order to reproduce the location and the 
amplitude of the maximum of the transmitted pulse in the celebrated experiment by Hau 
et al. [7]. We then get aout/(^in ~ 1-6, a broadening consistent with the observations, and 
^out ~ 0.16. The asymmetry being very slight, the delay of the maximum is very close to 
Tg and Coutit) is well fitted by the normal (Gaussian) form el^)-{t). A perfect fit is obtained 
by using the improved form e),J^f(t). The contribution eoffit) associated with hoff{t) is 
actually too small to be visible. Conversely hoff{t) will be responsible for the generation of 
short transients when ej„(t) comprises localised defects. As expected and recently discussed 
about the EIT experiments j27|, the front of these transients will propagate at the velocity 
c (instantaneously in our local time picture). Their peak amplitude will be especially large 
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FIG. 8: Propagation of a Gaussian-shaped pulse in the EIT experiment. The parameters are as in 
Figdwith cjj„, = 1.5 fxs. The fuh and dashed hues respectively are the exact intensity-profile of the 
pulse and the normal form. The location of the two maximums differ by 0.2 fj,s in agreement with 
the relation Atout = o./^'^out- Inset : the corresponding pulse-envelopes for —5 fis < t < 15 fis . 




when the defects consist in discontinuities. Consider again a square-shaped incident-pulse, 
the total duration 2Tp and the amplitude rj of which are such that its area and its variance 
equal those of the Gaussian-shaped pulse (FiglHI). eo//(t) is then easily derived from ho f fit). 
It reads 

eo//(t)=r/[/(t + r,)-/(t-r,)] (13) 

with 

f{t') = U{t') - 7,.A 1^ ^'^j'l^^^'f exp i-lbax) dx. (14) 

Each discontinuity in ej„(t) actually originates a large transient. Its initial amplitude is 
equal to that of the incident pulse and its successive maximums of intensity occur at the 
instants jf^/2A'yi,a later, jin being the n^^ zero of Ji{u) Note that the amplitude of the 
transients exceeds that of the smooth part of eout{t)- In a real experiment however the finite 
values of the rise and fall times of the incident pulse and of the detection bandwidth will 
generally limit the importance of the transients. By a deliberate reduction of the detection 
bandwidth, it is even possible to bring their intensity to a very low level without significantly 
affecting the delayed Gaussian-like part (see inset of FigE)). 

The results obtained on the model EIT-arrangement hold for an extended class of sys- 
tems having a transparency-window in a wide absorption profile. They only lie on three 
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FIG. 9: Propagation of a square-shaped pulse in the EIT experiment. Parameters as in the two 
previous figures. The upper and lower curves are the envelopes of the incident and transmitted 
pulses. The peak amplitude of the transients is slightly smaller than its theoretical value due to 
the finite time-resolution of the computations (0.4 ns). Inset : the first transient expanded on a 
0.2fj,s time-interval (bottom) and the same (xlO) after passage through a f^* order filter of time- 
constant o"out/100 (top). Such a filtering does not significantly affect the smooth part of the pulse 
which keeps very close to that obtained with an incident Gaussian-shaped pulse but reduces the 
amplitude (resp. the intensity) of the transient by a factor of about 30 (resp. 900). 



assumptions: (i) is Lorentzian in the far wings of the absorption profile (ii) the opaque 
regions are much wider that the transparency-window (iii) the transfer-function does not 
significantly deviate from the normal form in the transparency-window. The first condition 
(i) is generally met even when T{Q) is not Lorentzian in its central part. Anyway, it is not 
essential. If it is not met, the detailed shape of the transients is modified but not their 
main features (instantaneous transmission, duration proportional to the inverse of the spec- 
tral width of the opaque regions). The conditions (ii) and (iii), which are closely related, 
are met in the EIT experiments when the medium transmission is good at the frequency 
of the optical carrier (see before) but this is not always sufficient. As a counter-example 
we consider the experiment achieved by Tanaka et al. in an atomic vapour with a natural 
transparency- window between two strong absorption lines [2^. The complex gain-factor 
reads 

1 11 



A 
"2 
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where 2A is the doublet sphtting. Despite an apparent similarity, the associated transfer- 
function dramatically differs from that obtained in EIT when Qs > {iba — lea)- Indeed 
the two involved lines are here purely Lorentzian (not hybrid) and a good transparency 
at n = is achieved only if A ^ 7. We then get Hq ^ exp (— A7^/A^), Tg ^ y47/A^ , 
cr^ ^ QA'y'^/A'^, a ~ — 6y47/A^ and ^ ~ — A^/(7^v^6]4). Choosing the physical parameters 
such that Hq, Tg and cr^ equal their values in the EIT experiment , we actually obtain a 
quite different gain profile with opaque regions whose width is smaller than that of the 
transparency- window fFigirUj). This entails that the approximation ^ \H^'^\il)\ 
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FIG. 10: . Gain-profile |i?(r2)| for the two absorption-lines arrangement (full line) compared to 
the corresponding normal gain-profile |-fr^^(r2)| (dotted line). Hq, Tg, a and frequency unit as in 
FiglTJ Inset : envelope of the transmitted pulse for a Gaussian-shaped incident pulse (full line) 
compared to that obtained with the EIT-arrangement (dotted line). o"i„ = 1-5 [is 



only works in the immediate vicinity of i7 = 0. It is the same for the law $(^^) oc il, the 
skewness being very large ~ —7.6 ). In such conditions, even a Gaussian-shaped pulse is 
strongly distorted (see inset of FiglTUI). In fact the two absorption-lines arrangement allows 
one to attain large fractional delays Tgjoin with moderate distortion and broadening but this 
requires to involve much larger absorption-parameters and to accept a lower transmission. 
For example, Tanaka et al. succeeded in obtaining Tg/ain ~ 13 with Gaussian-like pulses (see 



Fig. 4c in |29|) but the peak-intensity of the transmitted pulse was 75 times smaller than that 
of the incident pulse. Their results are well reproduced with our two-lines model by taking 
A = 2.6 X 10^ and A/7 = 110. We then get a'^/af^ ~ 1/25 ^ 1 and we actually are in a case 
of low distortion as previously discussed (see FigEJ. However, pulses with discontinuities 
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are excluded. Since H{oo) ^ Hq, the resulting transients (not delayed) would indeed be 
much larger than the delayed part of the transmitted pulse and would obscure it. Moreover, 
due to the narrowness of the opaque regions, the time scales of the transients and of the 
delayed part do not considerably differ and it is thus impossible to filter out the former 
without denaturing the latter. 

VI. SUMMARY AND DISCUSSION 

Privileging the time-domain analysis, we have studied the linear propagation of light 
pulses, the frequency of which coincides with that of a pronounced maximum in the trans- 
mission of the medium. An important point is that substantial pulse-delays are only attained 
when the corresponding transmission exceeds the minimum one by a very large factor C 
(contrast). The impulse response of the medium then tends to a normal (Gaussian) form, 
irrespective of the line profile associated with the transmission peak. The propagation of 
arbitrarily shaped light-pulses with significant delays and low distortion is possible when the 
rms duration cr of the medium impulse-response is small compared to that of the incident 
pulse ((Tin), which should itself be small compared to the group delay Tg. The fulfilment 
of this double condition requires systems where C is extremely large, typically several tens 
of thousands of dB in a logarithmic scale. Systems with such contrast have actually been 
used j2| but, in most slow-light experiments, C ranges from 30 to 600 dB. Significant frac- 
tional time-delays Tgja^n keep attainable with such values of C by using incident pulses, the 
duration of which is comparable to or smaller than cr. As the medium impulse-response, 
the transmitted pulse then tends to acquire a normal (Gaussian) shape whatever its initial 
shape is. This reshaping is particularly striking when the incident pulse is square-shaped 
but is reduced to a simple broadening when the latter is itself Gaussian-shaped. Despite 
its asymptotic character, the normal form generally provides a good approximation of the 
shape of the transmitted pulse. More precise shapes are obtained by a perturbation method, 
allowing us in particular to specify how much the delay of the pulse-maximum deviates from 
the group delay. 

All these results have been first established by assuming that the transmission peak is 
put on an uniform background. We have shown that they also apply when the transmission 
peak is associated with a transparency window in an absorption-profile of finite width. 
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This however requires that the nearly opaque regions flanking the transparency window be 
considerably wider than the latter. Other things being equal, there are then no differences 
between the cases of uniform and non uniform transmission-backgrounds, at least when the 
envelope of the incident pulse is smooth. Conversely localised defects in this envelope will 
be responsible for the generation of very short transients which complement the normal 
(Gaussian) part of the signal. The front of the transients is instantaneously propagated in 
our local time picture (that is at the velocity c in a dilute sample). In extreme cases, their 
amplitude may be comparable to that of the delayed signal but, due to their location and 
their duration, they can easily be eliminated without altering the latter. 

The slow and fast light experiments have a common feature. In both cases, the observation 
of significant effects requires media with a very large contrast between the maximum and the 
minimum of transmission. This results from the causality principle and implies severe limits 
to the effects attainable in fast-light experiments, whatever the involved system is 3^ . From 
this viewpoint the slow-light case is obviously less pathologic and the constraints, although 
real, are much softer. 
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